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1.0 SUMMARY 


This report presents the results of a combined experimental (wind tunnel test 
results) and theoretical analysis utilizing random harmonic analysis 
techniques to predict the dynamic response and the structural dynamic loads of 
flat plate photovoltaic arrays due to wind turbu(ence. Guidelines for use In 
predicting the turbulent portion of the wind loading on future similar arrays 
using the results of this study are presented. 

The dynamic response and the loads dyn^Juic magnification factor of the two 
array configurations (a four post array and a two post array) are similar. 
Figures 1-1 and 1-2 are typical magnification factors for the two array 
configurations located In Identical positions within the array field. The 
figures show the magnification factors at a mid chord and outer chord location 
on uhe array Illustrated and at four points on the chord. The wind tunnel 
test experimental rms pressure coefficient that the magnification factors are 
based on Is also shown on each figure. In general, the largest response and 
dynamic magnification factor occur at a mid chord location on an array and 
near the trailing edge. A technique employing these magnification factors and 
the wind tunnel test rms fluctuating pressure coefficients to calculate design 
pressure loads due to wind turbulence Is presented In Figure 1-3. 
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Experimental RMS Pressure Coefficient 




< r 

f-igure 1-1. Configuration I - Dynamic Analysis Results for Sth Array 
with an Array Field 
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1st Mode Natural Frequency Hz 1st Mode Natural Frequency 



Experimental RMS Pressure Coefficient 



Figure 1-2. Configuration // • Dynamic Analysis Results for 5th Array 
within an Array Field 
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1st Mode Natural Frequency ^ Hz 1st Mode Natural Frequency 



step A (Decisions) 


( 1 ) 


Design configuration 
matches ? 

Configuration I (See 

II Figure 3-1) 


( 2 ) 


Design configuration 1st 
structural natural 
frequency excited by wind 


(3) 


Rms unsteady pressure 
coefficients for 
appropriate design 
configuration - Appendix 8 


(4) 


Level of probability not 
to exceed 
68.27% = 1 o 
95.45% = 20 
99.73% = 30 


(5) 


Design wind speed at 
an aleV'ition of 10 meters 



Obtain magnification factor 
(M.F.) from appropriate 
figures 3-3 to 3-8 using 
decisions 1 and 2 



Multiply M.F. by rms 
unsteady pressure 
coefficient from decision 3 
* design pressure 
coefficient with probability! 
not to exceed 68.27% 



Multiply Step C by level 
of probability not to 
exceed from decision 4. 




Multiply Step D by dynamic 
wind pressure from decision 
5 * design pressure load on 
array surface due to 
non-steady portion of wind 


Figure * 3. Flow Chart to Calculate Array Design HllirHl Loads due to 
the Non-Steady Portion of the Wind 
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2.0 INTRODUCTION 


This report summarizes a combined theoretical experimental analysis of the 
structural dynamic loading on long, flat plate photovoltaic arrays resulting 
from exposure to the non- steady portion of the wind environment. This report 
is an extension to the analyses reported In OOE/JPL 954B33-79/2 and -81/3. The 
study was performed under contract number 954833 to the Jet Propulsion 
Laboratory as part of the Engineering Area Task of the Low Cost Solar Array 
(LSA) Project. This project Is being managed by JPL for the Department of 
Energy, Division of Solar Technology. 

2.1 Objectives 

The cost associated with the design and construction of solar photovoltaic 
arrays to produce electric energy from sunlight Is an Important factor In the 
acceptance and use of solar energy. The load due to wind on an array and on 
Its support structure strongly Influences the design and ultimately the cost of 
the photovoltaic panels, panel and array support structure and foundation of 
the array. It Is, therefore, essential to determine the true maximum wind load 
tnat the array will experience during Its lifetime In order to minimize the 
structure costs. The objective of this study was to establish wind load guide- 
lines on flat plate photovoltaic arrays for that portion of the wind 
environment that Is considered non- steady (turbulence). 

2.2 Discussion and Background 

Three factors affect the amount of wind loading on a body: the fl field in 

which the body Is placed, the aerodynamic characteristics of the body Itrelf, 
and the dynamic response of the body due to the wind loading. Altl ugh the 
structural loads resulting from this latter factor are not totally composed of 
aerodynamic forces (they also Include Inertia forces), these structural loads 
do result from the wind loading, or more precisely, the fluctuations In wind 
loading. 

The flow field of the type that would be found around arrays In an array field 
has three aspects: 1) the steady state flow before It encounters any 

obstacles, 2) atmospheric gusts, and 3) turbulence. The steady state flow 
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is characterized by of a shear layer adjacent to the ground whose shearing 
effect?; decrease with elevation above the ground until a unlfortn flow is 
attained. Gusting is the result of velocity variations and changes in the 
direction of the prevailing wind due to atmo^phei ic instabilities. Turbulence 
may be caused by several factors. Gusting can cause turbulence when adjacent 
volumes of air are moving at different velocities, thus producing a shearing 
effect. Terrain roughness causes turbulence because of shearing effects. An 
obstacle in the path of the flow can also create turbulence by upsetting the 
flow and causing eddies and vortices to form. In addition, the shape of the 
body will affect the characteristics of the turbulence. Turbulent flow is 
highly complex, with varying frequencies and intensities occurring .n a random 
manner. 

The effect on the array forces due to the flow field is a function of the 
airflow characteristics and the resulting pressure distribution over the array. 
When the flow is turbulent, the pressure distribution and, consequently, the 
forces exerted on the array are nonuniform in frequency and intensity. These 
forces may cause vibrations in the structure resulting in additional structural 
dynamic forces on the array. Since turbulence varies in frequency and 
intensity, the resulting loading will also vary as a function of the frequency 
and intensity. 

2.3 Study Requirements 

The requirements of this study involve analysis and test. They are: 

1. Wind tunnel test and data reduction in the form of auto and cross 
spectrums. 

2. Theoretical structural dynamic modeling of arrays. 

3. Dynamic analysis results. 

4. Establishment of design guidelines for estimating unsteady wind loads 
on photovoltaic arrays. 
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The foll^lng Is a summary of the statement of work for Phase IV. 

1. Array Vlb.atlon Characteristics 

i) A structural finite element computer model will be developed in 
sufficient detail to adequately define the structure for usu in the 
calculation of structural mode shapes of a photovoltc:ic array. The 
tilt angle of the array will match one of the tilt angles of the 
arrays used in the Phase III study (reported in reference 1). 

ii) Structural vibration characteristics of the array will be calculated 
to produce mode shapes and ge ^alized mass and stiffness data for 
the array model developed in i). 

2. Array Structural Dynamic Load 

i) The generalized mass and stiffness data and the vibration mode shapes 
will be combined with the steady state aerodynamics from one array 
field configuration chosen from the Phase III test data to generate 
the structural dynamic equations of motion. 

ii) The Phase III wind tunnel pressure time history data will be analyzed 
for the array configuration chosen in 2 (i) to yield auto- and cross- 
spectral density functions and will be used in calculating the 
generalized forcing functions for tne equations of motion. 

3. Structural Dynamic Loads 

Structural dynamic responses and loads for the arrays will be 
calculated from the equations of motion and load equations developed 
in 2. The resulting aerodynamic forces from the dynamic analysis 
will be compared to the steady state aerodynamic forces to indicate 
the level of increase to be expected from th? dynamic forces. 

2.4 Report Organization 

The remainder of this report presents the results of the dynamic analysis and 
its formulation. Section 3.0 presents the basic technical approach, details 
the results and discusses the pertinent findings. Conclusions are presented in 
Section 4.0. The use of the results as design guidelines are presented for 
wind turbulence in Section 5.0. New Technology and References are outlined in 
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Sections 6.0 end 7.0 respectively. Appendix A presents the analysis procedure; 
Appendix B, the wind tunnel Measured ms delta pressure coefficients; and 
Appendix C» the Colorado State University test report for the test and test 
data reduction pertinent to the dyniMic analysis. 


t)uli;iNAL PAGE IS 
r\v P(y^H i)VAJ m’ 
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3.0 TECHNICAL APPROACH. RESULTS AND DISCUSSION 


Because the pressure fluctuation on the arratys caused by the Mind turbulence 
Is random In nature, random harmonic analysis techniques can be applied to 
determine the dynamic response of the arrays. The technical approach Mas to 
develop theoretical mathematical models of arrays and solve for the response 
of the arrays using Mind tunnel measured pressure fluctuations on the arrays 
as the forcing function. Pressure measurements mere recorded simultaneously 
at several locations on rigid arrays during a Mind tunnel test conducted at 
Colorado State University. These pressure data Mere then reduced to yield 
auto and cross spectra at and betMeen each pressure tap. The test procedures 
and results of the wind tunnel test Is presented In Appendix C. 

The structural dynmnic response of an array Is dependent on the structural 
characteristics of the array as well as the forcing function. For this study, 
tMO typical arrays (Figure 3-1) Mere modeled theoretically utilizing the modal 
approach; the mode shapes are shOMn In Appendix A. The first configuration 
Mas a four post array designed by JPL (see Ref. 5 for the detailed design 
configuration) and the second configuration Mas a tMO post array (support 
posts at the center of each end of the array). The theoretical equations for 
the dynamic response of the arrays Mere calculated and Included the equation 
of motion and pressure load equations on the arrays. The theoretical 
procedure Is detailed In Appendix A. The fluctuating Mind Induced pressure 
forcing functions (auto and cross spectra*) on these arrays Mere simulated 
from the Mind tunnel test for the wind conditions on the first array of an 
array field exposed to the vrind, the first array protected by a fence, and the 
fifth array Mithin an array field. The Mind Mas considered to be head-on to 
the arrays. 

The Mind tunnel test results consisting of auto and cross spectra of the array 
pressure coefficients are of significant Importance to the dynamic response of 
the arrays. Figure 3-2 shOMS a typical nonnallzed cross spectrum for the 5th 
array located in an array field Mith a tilt angle of 35® and a 90 mph Mind. 

Of importance In the cross spectrum Is that the phase angle in all cases Mere 


*The cross- spectrum relates pressure versus frequency at two different points 
(1 and j); the auto- spectrum Is similar except at a single point (1=j). These 
terms are used In the analysis as shOMn In Appendix A. A detailed technical 
de<c>’ip( ion of the auto and cross spectrum can be found In References 3 and 4. 
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Configuration I 



Height ■ 977 lbs 

First Plate Bending Mode ■ 10 Hz 


Configuration 11 


|d 12 ft H 



Weight • 933 lbs 
Pitch Mode - 3 Hz 
First Plate Bending Mode - 10 Hz 


FigunS'l. Array Configuntion Chtncttrhtks 
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zero. That 1s» the pressure at each pressure tap Is in phase with each other 
pressure tap on the array. The In^rtance of this feature Is that only 
synmetric Modes of vibration about the center line of the array can be excited 
by the Mind turbulence. Another important aspect of the spectra is that for a 
ninety mph wind the power in the wind as indicated by the array pressure 
coefficient spectra becomes insignificant (<1X of maxiMum) at frequencies 
greater than 5 Hz. 



r equency ^ Hz 


Figun 3-2. Typkst Cross Spsetrum bstwstn hwssun Tsps on sn Array 
in a 90 mph Wktd 


Utilizing randotn harmonic analysis techniques and the properties of auto and 
cross spectrums, a structural dynamic analysis was performed for the two array 
configurations. The results presented as magnification factors defined as the 
theoretical RMS dynamic load divided by the experimental RMS pressure load is 
shown in figures 3-3 to 3-5 for configuration I and figures 3-6 to 3-8 for 
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configuration II. Figures 3-3 to 3-5 show the effect of location of 
configurate I arrays In an array field (1st arr^ without a ^ence, 1st array 
protected by a fence, and 5th array within an array field, respectively). 

Each figure also shows the effect of the first modal structural frequency on 
the magnification factor and the effect of location on the array (mid chord 
.r .d outer chord locations). In general, the largest magnification factors are 
ocated near the trailing edge of the array (downwind edge) and at the mid 
chord location. The magnification factor decreases with an increase in the 
virst mode structural frequency. For the condition where the first mode 
frequency is between one and two cycles/sec, the magnification factors change 
^ ipidly. This is caused by the combination of the very flexible structures 
3S indicated by the low natural frequency) and the aerodynamics interacting 
to reduce the stability of the structure. At frequencies above 3 Hz, the 
structure becomes stiff enough so that the interaction with the aerodynamics 
does not present a problem. The magnification factors for the 5th array and 
the 1st array protected by a fence tend to be slightly larger than the 1st 
array that is unprotected. This is expected since the fence and the upwind 
arrays will increase the level of turbulence in the wind compared to the wind 
flowing over the 1st array. 

The results for configuration II arrays (figures 3-6 to 3-8) are similar but 
of reduced ■ jgnitude compared to the results for configuration I arrays. The 
reduced magnitude is a result of the larger weight per unit surface area tha". 
is required for static load considerations by configuraticr. II compared to 
cdfiguration I. As a result of the higher weight, less dynamic response 
occurs for the same 1st mode structural frequency. Also, the magnification 
factors for the outer chord and the mid chord locations are very similar in 
shap<.! and mag itude for configuration II. 

In n r.eral, the structural configuration (weight, stiffness and shape) will 
af.act the dynamic response of the structure. The larger the weight and 
stiffness, the less the dynamic response for a given forcing function. The 
weiv-ht and stiffness of the ar»'ay considered as configuration I was about as 
light and rlexible as possible and yet maintain structural integrity for the 
rtc*'i' loads. The naninal first structural natural frequency of this array 
that would be excited by the wind was just slight below 10 Hz. For both 
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configurations, the structure whose 1st wind excited structural frequency Is 
below 5 Hz Is sufficiently flexible that It Is questionable whether It could 
withstand the static loads with sufficient confidence to be a viable 
structure. Because for configuration I to have a lowest natural frequency of 
5 Hz requires the structure to be approximately 4 times as heavy or 4 times as 
flexible as the nominal structural configuration, the results for the 
configuration with a 5 Hz frequency should be very conservative. If 5 Hz Is 
considered as the lower limit for a variable structure, the maxiroimi 
magnification factor for configuration I Is 1.3 and configuration II Is 1.08. 
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Experimental RMS Pressure Coefficient 



f ijw/e 3-3. Configuration / • Dynamic Analysis Results for 1st Array 
Unprotected from the Wind 
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1st Mode Natural Frequency^ Hz 1st Mode Natural Frequency Hz 




Experimental RMS Pressure Coefficient 



Figun 3-4. Configuration I - Dynamic Anaiytis Ratuits for 1st Array 
Protactad from tha Wirtd by a FarKa 


15 


1st Mode Natural Frequency Hz 1st •‘•tode Natural Frequency 




Experimental RMS Pressure Coefficient 



Figure 3-5. Configuration / • Dynamic Anaiysis Results for 5th Array 
within an Array Field 
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1st Mode Natural Frequency hz 1st Mode Natural Frequency 




Experimental RMS Pressure Coefficient 





Figure 3-d ConfigunMn II ■ Dynamic Analyw ffasults for 1st Array 
Unprotactad from ihe Wind 
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1st Mode Natural Frequency 1st Mode Natural Frequency 



Experimental RMS Pressure Coefficient 



Figun 3‘7. Configuration II • Dynamic An^ysis Results for 1st Array 
Protected from the VUInd by a Fence 

Ifl 


1st Mode Natural Frequency ^ Hz 1st Mode Natural Frequency 



Experimental RMS Pressure Coefficient 



a 

u 


< x: 

Figun 3-8. Configuration II - Dynamic Analysis Results for 5th Array 
within an Array Field 
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1st Mode Natural Frequency ^ Hz 1st Mode Natural Frequency 


4.0 CONCLUSIONS 


Although the structural configuration affects the dyncunic response due to wind 
turbulence, the general trends of the two configurations studied are similar. 
The mid chord location on an array has the largest response and the largest 
magnification factors on be h configurations. The outer chord magnification 
factor of configuration I Is considerably lower than the mid chord factor 
because of the corner constraints, whereas for configuration II, the outer 
chord magnification factor Is only slightly reduced compared to the mid chord. 
The largest dynamic response occurs near the downwind edge of the array and 
decreases as the 1st structural modal frequency Increases. to minimize 

the dynamic loads due to turbulence requires the maximum pr j1 stiffness 
for the structure. A dynamic magnification factor of 1.3 ap to be a 
conservative factor for any type of flat plate photovoltaic array 3truc^ure. 

It should be noted that for the types of structures used In this analysis, if 
the 1st wind excited structural frequency Is below 5 Hz the structure Is 
sufficiently flexible that It Is questionable whether It could withstand the 
static loads with sufficient confidence to be a viable structure. For 
configuration I to have a lowest natural frequency of 5 Hz requi*‘as the 
structure to be approximately 4 times as heavy or 4 times as flexible as the 
nominal structural configuration, thus, the results for the configuration with 
a 5 Hz frequency should bc: very conservative. If 5 Hz Is considered as the 
lower limit for a variable structure, the maximum magnification factor for 
configuration I Is 1.3 and configuration II Is 1.08. 


MM lIlIENIMijjiUd: 
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5.0 ARRAY DESIGN GUIDELINES FOR WIND TURBULENCE 


Based on the results of the boundary layer wind tunnel test and the results of 
this dynamic analysis, the following design guidelines are given for 
determining wind loading on photovoltaic flat plate arrays for the turbulent 
portion of the wind. These guidelines are valid fo” arrays at least 2 chord 
lengths from the edge and for a wind velocity profile that approximates a 1/7 
power law. However, because of the edge constraints of the arrays In a 
dynamic sense, and because of the conservatism In this analysis, the wind 
designs loads calculated using these guidelines would probably be also 
adequate for the side edge arrays. 

FI !re 5-1 presents guidelines to calculate design pressure loads on array 
surfaces due to the non-steady po**t1on ot the wind. (The deslpn guidelines 
for the steady state portion of the wind that needs to be combined with the 
turbulent po» cion for the total design load was documented In reference 1.) 
This procedure requires a number of decisions which Includes; 1) the choice of 
confl duration I or II must be representative of the design under 
consideration; 2) the choice of the rms unsteady pressure coefficients for the 
appropriate design array tilt cingle and array location from the wind tunnel 
results presented In Appendix B; 3) the level of probability not to exceed a 
given level (I.e., lo = 68.27%, 2a = 95.45%, 3a ^ 99.73%); 4) the design 
configuration first structural natural frequency that can be excited by the 
wind turbulence; 5) the design free stream wind speed. The appropriate 
magnification fact from figure 3-3 to 3-8 multiplied by the rms unsteady 
pressure coefficient chosen from Appendix B is the design rms pressure 
coefficient for wind turbulence with a 68.3% probability of not exceeding. 

For a higher probability such as 95.5%, this design rms pressure coefficient 
would be multiplied by 2.0. The following is an example oopiying this 
procedure. 

Assume that the design Is similar to configuration I, has a tilt angle of 35“ 
and the first structural fn ,uency normal to the array is 5 Hz. For any 
interior arrays, the rms unsteady pressure coefficient for this condition 
extracted from Appendix B is roshown in figure 5-2. The magnification factor 
(M.F.) for 5 Hz at the mid span location (figure 3-5) and these M.F. 
multiplied by the rms ACp from figure 5-1 is shown in table 6.1 as a function 
of chord position. Assuming that a level of probability not to exceed 95.5% 
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Step A (Decisions) 


Design configuration 
matches ? 

Configuration I (See 

II Figure 3-1) 


( 2 ) 


Design configu<^ation 1st 
structural natural 
frequency excited by wind 


(3) 


Rms unsteady pressure 
coefficients for 
appropriate design 
co nfiguration - Appendix B 

(4) 


Level of probability not 
to exceed 
68.27% = lo 
95.45% = Zn 
99.73% = 30 


(5) 


Design wind speed at 
an elevation of 10 meters 



Obtain magr.ification factor 
(M.F.) from appropriate 
figures 3-3 to 3-8 using 
decisions 1 and 2 



Multiply M.F. by rms 
unsteady pressure 
coefficient from decision 3 
= design pressure 
coefficient with probability 
not to exceed 68.27^; 



Multiply Step C by level 
of probability not to 
exceed from decision 4. 



Step E 


Multiply Step D by dynamic 
wind pressure from decision 
5 = design pressure load on 
array surface due to 
non- steady portion of wind 


Figure 5-1. F.’ow Chart *o Ca/cuia*f. Array ri- r,rdueto 

the Non-Steady Portio.i 
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of the design load due to turbulence is required, the rms ACp times the M.F. 
is increased by a factor of 2. The design delta pressure coefficient for the 
turbulent portion of the wind is shown in table 5-1 and in figure 5-3. The 
pressure due to a 90 r.ph wind is also shown in table 5-1 and figure 5-3. 

The design load for the front array (either protected or unprotected by a 
fence) is calculated in the same manner using the appropriate rms delta 
pressure coefficient and magnification factors. This procedure can be 
simpified by utilizing the highest M.F. (M.F. » 1.3 for modal frequencies > 5 
Hz) for all conditions. This would eliminate the need to determine the design 
configuration natural frequency as well as not requiring the determination of 
the M.F. along the chord from figures 3-3 to 3-8. 
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.2 .4 .6 .8 1.0 

Fraction of Chord (S/C) 



Figun S-2. Wind Tunnni Tm RMS Pmsun Coefficients on dteSth Amy 
in en Amv Field 
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6.15 > 
P 

(psf) 

4.1 • 


2.05 ^ 



Figure O'S. Exempie Design Wind Loeds on Interior Arreys in an Array 
Field Using the Design Guidelines 
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TtbhS-l. RmulUforan IntariorArmy Uting ttm Dmign GmMiim 


Steps 

see 

Figure 5-1 


S/C 

- 

.96(T.E) 

.74 

.26 

.04(L.E: 

Step A 

Rras ACp (see Figure 5-2) 

.104 

.095 

.133 

.189 

Step A 

q (psf) (90 mph wind) 

20.5 

— 

— 

— 

Step B 

M.F. (see Figure 3-5) 

1.19 

1.27 

1.17 

1.07 

Step C 

(M.F, ) ( rms ACp) 

.1238 

.1207 

1 

.1556 

.2022 

Step D 

Design ACp 

(2o ) (M.F.) (rmsACp) 

.2475 

.2413 

.3112 

1 

.4045 

Step E 

Design Pressures 
(psf) 

5.07 1 

i 

1 

4.95 

6.38 

8.29 


T.E. = Trailing Edge 
L.l - Leading Edge 
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6.0 NEW TECHNOLOGY 


No reportable Items of new technology have been Identified by Boeing during 
the contract of this work. 
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Appendix A - Analysis Procedure 


A mathematical model of the photovoltaic array included the important 
structural, mass and aerodynamic characteristics necessary to simulate the 
dynamic response of the array to the turbulence generated by the upwind arrays 
and that contained in the free stream wind. 

The structural model included the elastic characteristics of the complete 
array. A finite element program NASTRAN was used to model the structure and 
to calculate modes of vibration. The mode shapes for the two configurations 
that can be excited by aerodynamic turbulence are shown in Figure A-1. 
Generalized aerodynamic forces representing the aerodynamics on the array were 
computed using the steady state pressure distributions measured in Phase III 
(reference 1 ) and combined with the modes of vibration to produce the 
following equations of motion for use in a Boeing developed Ranuom Harmonics 
analysis computer program (reference 2 ). 

[Mi]|q|+ [Majjqj + [ M 4 ] | q + [ M 5 ] | q j = [c]|apJ 

generalized structural stiffness and structural damping matrix 
generalized mass matrix 
aerodynamic stiffness matrix 
aerodynamic damping matrix 
vector of generalized coordinates 

matrix of generalized forces associated with unit pressures at each 
of the excitation points 
vector of excitation pressures 

A corresponding set of load equations using the force summation method was 
developed in the form 

|l . [»i]|‘)j t[M3]|qj+[«4]|qj<^[«5]|5j+[c]|Ap| 

Where is the load vector and the [ Mi ] matrices are as defined for the 
response equations, except that they represent generalized load coefficients. 
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The auto and cross power spectral densities of the pressures measured in the 
wind tunnel was used to excite the array mathematical model. The loads output 
spectrum was calculated utilizing the following equation: 

i=l j=l 

where 

X-i (oO transfer function for system response to excitation i 
• denotes conjugate 

^ij i = j auto power spectral density for pressure i 

i ft J cross power spectral density between pressures i and j 

Integration of the output spectrum produces the rms load response 

_ / x1/2 

L = ( / L(u))do) ) 

' 0 ' 

Because the pressure taps on the wind tunnel model in the Phase III study were 
located on only one spanwise station, other spanwise stations were assumed to 
have an identical auto- and cross- spectrum chordwise and a correlation 
function of unity spanwise. The assumption of the same chordwise spectrum at 
each spanwise station is a valid assumption since the process is random. 
However, the assumption of a correlation function of unity between spanwise 
pressure taps is conservative. 
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Mode 1 




Frequency Ratio: Mode 2/Mode 1 = 1.6 


FigunA’t. Array Mode Shapes 
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Mode 3 



Frequency Ratio: Mode 3/Mode 1 = 3.3 



FigunA'1. Amy Modt Sh9pes (Continumil 
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Configuration II 


I 


Frequency Ratio; Mode 2/Mode 1 » 3.6 



Frequency Ratio: Mode 3/ Mode 1 » 5.0 
FigunA-1. Amy Mod» Shapes (Concluded) 
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Wind Tunnel Test RMS Unsteady Delta Pressure Coefficients 

RMS pressure coefficients shown in Figures B-1 to B-3 of the fluctuating 
pressures were obtained from the wind tunnel test results reported in 
reference 1. These coefficients were obtained for three tilt angles (20”, 35" 
and 60"), and for the first array of an array field unprotected and protected 
by a fence, the fifth array within an array field, and for wind directions 
from the front and rear. The pressure coefficients were calculated based on a 
wind reference velocity at 10 meters and a 1/7 power law wind velocity 
profile. 
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Fraction of Chord (S/C) 


f igure B- 1. RMS Pressure Coefficients on the First Array 
of an Unprotected Array Field 
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PREFACE 


Simultaneous pressure measurements on flat plate photovoltaic solar 
arrays were conducted in a simulated atmospheric boundary layer charac- 
terized by a l/7th power law mean velocity distribution. This report 
describes random properties of fluctuating local pressures on solar 
arrays by cross-spectral analysis. The random properties to be analyzed 
include the auto- and cross-spectra for fluctuating pressures for several 
typical array configurations. 

The essential experimental configurations, including facilities, 
wind-tunnel models, instrumentation and the flow simulation technique, 
have been described in a preceding report [5] "Wind Pressures and Forces 
on Flat-Plate Photovoltaic Solar Arrays," Colorado State University 
Report CER80-81NH-JAP-MP-JEC13. The simultaneous pressure measurements 
presented in this supplementary report required some additional arrange- 
ments for instrumentation and data acquisition. These arrangements will 
be described in Sections 1 and 2 of this report. 



1.0 PRESSURE TAPS AND ARRAY CONFIGURATIONS 


1.1 Pressure Taps 

‘Hie simultaneous pressure measurements were obtained at four 
pressure taps along the chord on the upstream surface and four pressure 
taps on the downstream surface of the solar array. These pressure taps 
tested were taps 1, 4, 7 and 10 on the upstream surface, and taps 11, 

14, 17 and 20 on the downstream surface. The location of each pressure 
tap is shown in Figure 1 which duplicates Figure 10 of the preceding 
report. 

1.2 Array Configurations 

All array configurations tested for this phase of the study are 
listed in Table 1. Array locations are shown in Figure 2 which dupli- 
cates Figure 11 in the preceding report. Several typical array configu- 
rations are chosen to be presented and discussed in this report. For 
this study, the ground clearance and spacing of solar arrays were 0.25 c 
and 2.0 c respectively, and were not varied. The wind direction was 0” 
for all array configurations. The definitions of these test parameters 
are seen in Figure 2.. 

2.0 DATA ACQUISITION 

The mean wind velocity outside the simulated turbulent boundary 
layer in the Meteorological Wind Tunnel, shown in Figure 1 of the main 
report, was approximately 50 fps as was used for the preceding tests 
(giving a U^^^ z. 43 fps) . The outputs from the 8 pressure transducers 
were recorded simultaneously on digital magnetic tape for 35 seconds at 
500 samples oer second using a data acquisition system based on a 
Hewlett-Packard System 1000 minicomputer. The data werr then analyzed 
by the same computer. Pressure, normal force and pitching mome' t 
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coefficients presented in this report were referenced to the Been 
dynaaic pressure at 10 ■ elevation in the prototype. The power spectral 
densities of local peak pressures were evaluated with these pressure 
coefficients. 

2.1 Pressure and Noraal Force Coefficients 

Pressure and noraal force coefficients are respectively defined in 
Sections 3.2 and 3.3 of the preceding report. The sane definitions 
apply to this report: 


S 


P . AP 

^ref ^ *^ref 

c 

if ilCp(s)d« 

6 


( 1 ) 


2.2 Pitching Monent Coefficients 

Pitching moment coefficients, about the mid-chord of the solar 
array are defined by 




M 


( 2 ) 


Tef 


where M is the calculated pitching moment per unit surface area about 
the mid-chord of the array, is the reference dynaadc pressure of 

approaching flow at 10 m full-scale, and c is the chord length of the 
array. In order to calculate M, the eccentricity of the center of 
pressure, e, was obtained numerically by a curve fitting to the pressure 
distribution using linear interpolation schemes along the chord of the 
array . Thus 


M* ■ FN • e 


( 3 ) 


where FN is the normal force. 


^pitching moment lifting windward edge up is defined to be positive. 
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2.3 Spectral Analysis 

2.3.1 Matrix Arrangement 

A set of cross-power spectral densities of local pressure 
fluctuations on the solar array were calculated. For simplicity of 
notation, the pressure trp numbers 1,4,7,10,11,14,17 and 20 can be 
renamed as 1,2, 3 ,..., 8 in order. Any pressure tap number used in 
this report will be, hereafter, referred to as the renamed one. 

Matrices for the correlation and power spectral densities are 
written as follows. 


- 

1 

R. .(t) 
IJ 

II 

4. .(N) 
IJ 

II 


auto-correlation , 
cross-correlation, 

auto- spectrum , 

cross-spectrum , 


if 

if 

if 

if 


i = j 
i * j 
i = j 
i / j 


where t is time lag, N is frequency, and i,j refer to the pressure 
tap numbers. 


2.3.2 Correlation Functions 

For this analysis, auto- and cross-correlation functions are defined 


by 




l.i| 

-MO i 


1 

T-m» 


C (t) - c 
P4 p. 


■ Ifs. 

mean^ *'3 


(t+T) 


“^mean^ 


dx 


(4) 


where (t) and ®re instantaneous and time-averaged values 

of pressure coefficient with respect to tap number k. Defining a 


fluctuating component of pressure coefficient, C^(t), as 

* J\ 


C 

Pk 


(t) 




‘^mean. 


(5) 
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the auto-correlation function, at t « 0, becones for one tap 


R..(0) = C (t) C (t) 

11 Pi Pi 


= C* 


rms. 

1 


in which an overbar denotes a time averaging. 

2.3.3 Power Spectral Density Functions 
Power spectral density functions are defined by 


( 6 ) 


4. . (N) 
13 


0 


Rij(T) exp(-j2nNT)dT 


.2 


(7) 


where j appearing in the exponential function refers to j = -1, 
Tlie integral property of the auto-spectral function requires 


fJO 

1^. 


(N)dN = R..(0) = C 

11 11 p 


( 8 ) 


0 


rms. 

1 


Because of this property, the power spectral density can be normalized 
by 


«*.(N) = p- 

n C_ 


4 . . (N) 


with units 


P P 

*^nns . *^nns . 

1 1 


- [k] 


( 9 ) 


It is also conunon practice to normalize the frequency N by 
Nc 


N* 


U 


ref 


( 10 ) 


where 0^^^ is the reference wind velocity. 

The cross-spectral analysis was performed digitally by a Fourier 


Transform subroutine using standard techniques such as those in 
references [1] and [3]. Transforms were performed on 8 time segments, 
each 2048 samples in length, for each pressure record (recorded at 
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500 samples per second) . Transforms were combined to form cross-spectra 
and appropriate averaging across data points in frequency and time seg- 
ments was performed to reduce normalized standard error o' the spectrum. 
Because of memory limitations, the cross-spectra were only calculated to 
a maximum frequency of 125 Hz (N*— 1.0), This frequency retained virtual 1' 
all the energy in the fluctuating pressures. Normalized standard error 
of the cross-spectra reached a maximum of 11 percent at the lower 
frequencies. 

3.0 RESULTS AND DISCUSSIONS 

3.1 Pressure and Normal Force Coefficients 

Examples of time traces of normal force and pitching moment 
coefficients are shown in Figure 3 obtained with the first solar array 
at a = 35", without fence. The time-averaged mean and rms pressure 
coefficients are tabulated in Table 2 for six typical array configurations, 
including Runs 21331, 21344, 21346, 21352 and 21368. These values agree 
well with those found in the preceding report [5]. 

3.2 Power Spectral Analysis 

The auto-spectra lor each pressure tap are showr. in Figures 4 
through 9 for the six cases listed in Table 2. Comparing the auto- 
spectra for the individual pressure taps and for both surfaces of the 
solar array, the power spectra show that the energy is shifted to higher 
frequencies for the first array with fence and the fifth array, consistent 
with the smaller scales of turbulence expected behind the fence and within 
the array field. The rear side of the first array also shows this effect. 
The area unaer the auto-spectra over the available frequency range is 
compared to the measured variance in fluctuating pressure in Table 3. 
Figures 4 through 9 clearly show that the energy content in the power 
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spectra is concentrated at N* < 0.1. At N* 0.3, a spike in the 
power spectra is seen for both the upstream and downstream surf.xes. 

This spike in the power spectra is due to vibration of the wind-tunnel 
model. Moreover, for some cases, another spike at N* j:. 0.4 was obtained 
which is a second mode of vibration of the model. No attempt was made to 
model the full-scale structural stiffness or damping so that these 
response spikes do not indicate full-scale response. 

Both the real and imaginary parts of the auto- and the cross-spectra 
were calculated. Only the real part is shown in Figures 4 through 9. 

The imaginary nart was essentially zero for all cases. The implication 
of this finding is that the phase angles of the various frequencies are 
uncorrelated. This is a typical result in velocity or pressure 
fluctuations in turbulent bourdary layer flows. 

Figures 10 through 15 show cross- spectral plots for all combinations 
of taps for each of the six cases listed in Table 2. Cross-spectra for 
each array were grouped on plots so that similarly-shar ’ ■ ^rvos would 
appear on the same plot. The plots reveal th' same increase in energy 
at higher frequencies for arrays behind the fence and withir the array 
field as was observed for the auto-spectra. Also, the model natural 
frequency of vibration shows in the plots. 

4.0 CONCLUSIONS 

On the basis of the data presented in previous sections, the 
following conclusions can be drawn. 

1. Auto-spectra of local pressure fluctuations characteristically 
fall rapidly with increasing frequency. 

2 . Auto-spectra show higher energy at the higb»r frac^uencies 


where the pressure tap is within the array field or uck..fid a wind fence. 
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3. Cross-spectra showod properties similar to those in 1 and 2 
above. 

4. Cross-spectra between taps in s^arated zones were quite 
similar. 

5. The imaginary parts of both auto- and cross-spectra were 
essentially zero. 
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Figure 1. Position of Pressure Taps on Instrumented Model 
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Schematic Description of Array Field 
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Figure 3. Time Trace of the Normal Force and 
Pitching Moment Coefficients 
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Figure 5. Auto-Spectra for First Array at a = 60®, without Fence 
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Figure 7. Auto-Spectra for First Array at a ~ 160°, with Fence 
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Figure 9. Auto-Spectra for Fifth Array at a = 35®, without Fence 
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Figure 10. Cross-Spectra for First Array at a = 35", without Fence 
(continued) 
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Figure 11. Cross-Spectra for First Array at a * 60*, without Fence 
(continued) 
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Figure 12. Cross*Spectra for First Array at a * 145®, with Fence 
(continued) 
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Figure 12. Cross-Spectra for Fi 
(concluded) 
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Figure 13. Cross-Spectra for First Array at a = 160®, with Fence 
(continued) 
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Figure 14. Cross-Spectra for Fifth Array at a » 20®, without Fence 
(continued) 
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Table 1. List of Array Configurations Tested 


Run 

a 

Array 


Fence** 


Height 

(Hf/c) 

Distance 

(Xf/C) 

Porosity 

(%) 

21318 

120 

5 

— 

— 

— 

21321 

120 

1 

— 

— 

-- 

21323 

120 

1 

0.75 

2.5 

30 

21325 

145 

5 

— 

— 

— 

21327 

145 

2 

— 

— 

-- 

21329 

145 

1 

— 

— 

— 

21331 

145 

1 

0.75 

2.5 

30 

21338* 

145 

1 

— 

— 

-- 

21340 

160 

5 

-- 


— 

21342 

160 

1 

— 

— 

— 

21344 

160 

1 

0.75 

2.5 

30 

21346 

20 

5 

— 

— 

— 

21348 

20 

1 

-- 

— 

— 

21350 

20 

1 

0.75 

2.5 

30 

21352 

35 

5 

— 

— 

-- 

21354 

35 

2 

— 

-- 

— 

21356 

35 

1 

— 

— 

-- 

21358 

35 

1 

0.75 

2.5 

30 

21364* 

35 

1 

— 


-- 

21366 

60 

5 

-- 

-- 

-- 

21368 

60 

1 

-- 

-- 

— 

21370 

60 

1 

0.75 

..ft ft • 

2.5 

30 


*edge study (see Section 2.3 o^ the preceding report for the definition) 
**see Figure 2 for definition of parameters 
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Table 2a. Tiaie-Averaged Pressure Coefficients 


Run 

21331* 

Run 21344* 

Run 

21346 

a 

Array 

a Array 

a 

Array 

145 

1 

160 1 

20 

5 


Tap (original) 

Shnean 

^Pnns 

^Pmean 

^Prins 

Shnean 

^Prms 

1 (1) 

-.083 

.051 

-.105 

.052 

.020 

.051 

2 (4) 

-.068 

.057 

-.096 

.061 

.055 

.050 

3 (7) 

-.088 

.078 

-.133 

.072 

.075 

.069 

4 (10) 

-.107 

.068 

-.109 

.060 

.077 

.113 

5 (11) 

-.134 

.040 

-.096 

.040 

-.196 

.210 

6 (14) 

-.136 

.039 

-.096 

.038 

-.116 

.121 

7 (17) 

-.153 

.042 

-.121 

.042 

-.103 

.078 

8 (20) 

-.131 

.044 

-.107 

.045 

-.054 

.072 

‘with fence 
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Table 2b. Time-Averaged Pressure Coefficients 


Run 

21352 

Run 21356 

Run 

21368 

a 

Array 

a Array 

a 

Array 

35 

5 

35 1 


1 


Tap (original) 

^Pmean 

^Prms 

^Pmean 

Srms 

^Pmean 

^rms 

1 (1) 

-.013 

.060 

.087 

.054 

.223 

.089 

2 (4) 

.001 

.062 

.292 

.067 

.376 

.095 

3 (7) 

-.004 

• 

O 

00 

.436 

.107 

.482 

.148 

4 (10) 

-.002 

.120 

.451 

.178 

.279 

.177 

5 (11) 

-.140 

.147 

-.300 

.081 

-.417 

.066 

6 (14) 

-.104 

.109 

-.290 

.073 

-.398 

.063 

7 (17) 

-.126 

.089 

-.346 

.071 

-.472 

.063 


8 (20) -.092 .097 -.293 .072 -.417 .064 
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Table 3. Comparison of Integration of Auto-Spectra 


array = 

1, a = 35®, without fence 

pressure taps 

J 

^ 1 

(1. 1) 

0.990 

(2, 2) 

1.069 

(3, 3) 

1.16S 

(4, 4) 

1.C59 

(5, 5) 

1.004 

(6, 6) 

0.971 

(7. 7) 

0.986 

(8, 8) 

0.963 


array = 5, 

a = 35®, without fence 

(1, 1) 

1.035 

(2, 2) 

1.016 

(3, 5'r 

1.070 

(4, 4) 

0.986 

(5, 5) 

l.Obl 

(6. 6) 

o 

i 

(7, 7) 

0.991 

(8, 8) 

0.972 


*by theory, this quantity should be identically equal to 1 
(see Equation 8) 




